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Classical bifurcations and entanglement in smooth Hamiltonian system 



M. S. Santhanam^, V. B. Shcorey^ and Arul Lakshminarayan^* 
^Physical Research Laboratory, Navrangpura, Ahmedabad 380 009, India. 
Max Planck Institute for the Physics of Complex Systems, 
Nothnitzer Strasse 38., Dresden 01187, Germany 
(Dated: February 1, 2008) 

We study entanglement in two coupled quartic oscillators. It is shown that the entanglement, 
as measured by the von Neumann entropy, increases with the classical chaos parameter for generic 
chaotic eigenstates. We consider certain isolated periodic orbits whose bifurcation sequence affects a 
class of quantum eigenstates, called the channel localized states. For these states, the entanglement 
is a local minima in the vicinity of a pitchfork bifurcation but is a local maxima near a anti-pitchfork 
bifurcation. We place these results in the context of the close connections that may exist between 
entanglement measures and conventional measures of localization that have been much studied in 
quantum chaos and elsewhere. We also point to an interesting near-degeneracy that arises in the 
spectrum of reduced density matrices of certain states as an interplay of localization and symmetry. 

PACS numbers: 05.45.-a, 03.67.Mn, 05.45.Mt 



I. INTRODUCTION 

The study of entanglement is currently an active area 
of research in view of it being a physical resource for 
quantum information theory, quantum computing, quan- 
tum cryptography and telcportation [1]. At a classical 
level, entanglement does not have a corresponding coun- 
terpart. However, increasingly it is being realized that 
the nature of classical dynamics, whether it is regular or 
chaotic, affects entanglement in the quantized version of 
the system [2[. In general, larger chaos in the system 
leads to larger entanglement production. This has been 
established by considering kicked top models 3], bakers 
map [1], Dicke model [H, billiard in a magnetic field [6j], 
kicked Bose- Einstein condensates 0] and Ai"-atom Jaynes- 
Cummings model In contrast to these studies, the 
role of classical bifurcations in entanglement of chaotic 
systems has not received much attention. Even though 
entanglement is a purely quantum attribute, it is never- 
theless affected by the qualitative nature of the dynamics 
in phase space. The results to this effect are obtained 
primarily in the context of quantum phase transitions in 
the ground state of infinite systems in which the entan- 
glement is maximal at critical parameter values [l^ . For 
instance, for case of ions driven by laser fields and cou- 
pled to a heat bath, i.e, a form of Dicke model was shown 
to exhibit maximal entanglement of its ground state at 
the parameter value at which classical system bifurcates. 
Similar result for the ground state was reported from the 
study of coupled tops, a generalization of the two dimen- 
sional transverse field quantum Ising model |9j] as well 
from Jahn- Teller models . The ground state entangle- 
ment of mono-mode Dicke model is shown to be related 
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to Hopf bifurcation [Tl| . Qualitatively similar results for 
two component Bose-Einstein condensate are also known 
(Tst - In all these cases, the treatment is confined mostly 
to the ground state of the system that exhibits criticality 
and involves one single classical bifurcation. 



Do these results hold good for chaotic, smooth Hamil- 
tonian systems that do not exhibit criticality in the sense 
of phase transitions ? As opposed to a single bifurcation, 
what happens in bifurcation sequences where stability 
loss and stability gain interleave one another ? Both 
these questions explore the connection between chaos and 
entanglement in a physical setting that is different from 
the earlier studies. In the context of this work, we ex- 
amine a Hamiltonian system whose classical dynamics is 
controlled by a single tunable parameter. The changes 
in the parameter leads to changes in the phase space 
structure; for instance regularity to chaos transition and 
bifurcation sequences of fixed points. Typically, chaotic 
systems display a sequence of bifurcations. Consider, for 
instance, the coupled oscillator systems, a paradigm of 
chaos for smooth Hamiltonian systems and is related to 
atoms in strongmagnetic fields, the quadratic Zeeman ef- 
fect problems [ij]. In these cases, one particular sequence 
of bifurcation is a series of pitchfork and anti-pitchfork 
bifurcations p^ . The pitchfork corresponds to a peri- 
odic orbit losing stability and in the Poincare section this 
appears as a elliptic fixed point giving way to a hyper- 
bolic fixed point. The anti-pitchfork is when the periodic 
orbit gains stability. In this work, we consider coupled 
quartic oscillators and show that the entanglement in the 
highly excited states of the system is modulated by clas- 
sical bifurcations. We could place this in the context of 
works that lend support to the notion that for generic 
one-particle states there is a strong correlation between 
entanglement and measures of localization [l^, [13, [3 ■ 
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II. ENTANGLEMENT IN A BIPARTITE 
SYSTEM 

A pure quantum state l^") composed of many subsys- 
tems \(t)i) is said to be entangled if it cannot be written 
down as a direct product of states corresponding to each 
of the subsystem. 

I*)entangled ^ ® \h) ® \4>3) ®\4>n) (1) 

Thus, entanglement implies stronger than classical cor- 
relations. If p = IV') (■'/'I is the density matrix representa- 
tion for a pure state then the reduced density matrix 
(RDM) can be obtained by applying the trace operation 
to one of the degrees of freedom. Thus, 

pi = Tr2|^)(V'| p2 = Tri|7/^)(V'| (2) 

are two RDMs, whose one of the degrees of freedom is 
traced out. The notation Tr, denotes that the trace oper- 
ation is applied on the zth degree of freedom. Schmidt de- 
composition [l| provides a representation for ji/') in terms 
of product of basis states, 

N 

where \4>i)(i) and |0i)(2) are the eigenvectors of the RDMs 
pi and p2 respectively, and are the eigenvalues of either 
of the RDMs. The von Neumann or the entanglement 
entropy of pure state is given by, 

N 

S=-Y,K\og\, (4) 

1=1 

Thus, when S* = 0, the subsystems are not entangled and 
when > 0, they are entangled. The Schmidt decompo- 
sition provides a compact and unique representation for 
the given eigenstate (unique in the generic case when the 
non-zero spectrum of the RDM is nondegenerate) . 

III. HAMILTONIAN MODEL AND 
BIFURCATION SEQUENCE 

A. Quartic oscillator 

We consider the Hamiltonian system given by, 

H =pI+pI + x^ + y'^ + ax^y^ (5) 

with a being the tunable chaos parameter. For a = 
0,2,6, the system is classically integrable and becomes 
predominantly chaotic as a — > oo. This has been exten- 
sively studied as a model for classical and quantum chaos 
in smooth Hamiltonian systems [Toj and exhibits qualita- 
tively similar dynamics as the host of problems involving 
atoms in strong external fields. In the limit a — s- oo, it 



is also of relavance as model of classical Yang-Mills field 
(20| . To study the quantum analogue of this system, we 
quantize it in a symmetrized basis set given by, 

V'ni ,712 

{x, y) = 7V(ni , n2) [0„i {x)(j)n^ (y) + (j)^^ (x)(/)„i (y)] 

(6) 

where Af{ni,n2) is the normalization constant and 
(j){x)(f){y) is the eigenstate of unperturbed quartic oscil- 
lator with a = 0. The choice of this form of basis set is 
dictated by the fact that the quartic oscillator has 
point group symmetry, i.e., all the invariant transforma- 
tions of a square. Hence we have chosen the symmetry 
adapted basis sets as in Eq. [H] from Ai representation of 
Civ point group. 

Thus, the nth eigenstate is, 

{^,y) (7) 

j(ni,n2) = l 

where anj(ni,n2) 

= {tlj(x,y)\'ifn{x,y)) are the expansion 
coefficients in the unperturbed basis space. Note that 
ni,n2 are even integers and anj{rn,n2) = a„j(„2,«i) in 
Ai representation of Civ point group. The eigenvalue 
equation is solved numerically by setting up Hamiltonian 
matrices of order 12880 using 160 even one-dimensional 
basis states. 



B. Bifurcation sequence in quartic oscillator 

In a general chaotic system many bifurcation sequences 
are possible. However, a two dimensional Hamiltonian 
system can exhibit only five types of bifurcations p^ . 
One such prominent sequence is a series of pitchfork and 
anti-pitchfork bifurcation shown schematically in Fig [T] 
To reiterate, a pitchfork bifurcation takes place when a 
stable orbit loses stability and gives rise to two stable 
orbits. Anti-pitchfork bifurcations happen when a sta- 
ble orbit is spontaneously born due to the merger of two 
unstable orbits. We will focus on a particular periodic 
orbit, referred to as the channel orbit in the literature 
[2l[, given by the initial conditions {x, y — 0,Px,Py = 0}, 
which displays such a bifurcation sequence. The Poincare 
section in the vicinity of the channel orbit has interesting 
scaling properties and the orbit itself has profound influ- 
ence on a series of quantum eigenstates, called localized 
states, in the form of density enhancements or scars [2^ . 
Such density enhancements due to channel orbits have 
also been noted in atoms in strong magnetic fields or the 
diamagnetic Kepler problem ^] as well. 

The stability of the channel orbit in the quartic oscil- 
lator in Eq ([5]) is indicated by the trace of monodromy 
matrix J{a) obtained from linear stability analysis. It 
can be analytically obtained for the channel orbits [l^l 
as, 

Tr J(a) = 2\/2cos (^^1 + 4a) . (8) 

The channel orbit is stable as long as |TrJ(Q;)| < 2 and 
it undergoes bifurcations whenever TrJ(a) = ±2. From 
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Chaos parameter 

FIG. 1: The schematic of a typical bifurcation sequence 
involving a series of pitchfork (circles) and anti-pitchfork 
(square) bifurcations as a function of chaos parameter. The 
solid lines indicate that the orbit is stable and dashed line 
indicate instability. 




FIG. 2: The Poincare section for the quartic oscillator in Eq. 
[5] in shown for (a) a — 90 and (b) a — 90.5. Note that at 
a = 90 the periodic orbit undergoes a pitchfork bifurcation. 



this condition, it is clear that the bifurcations take place 
at Q!„ — n{n + 1), {n = 1,2,3....). Thus the channel or- 
bit undergoes an infinite sequence of pitchfork and anti- 
pitchfork bifurcations at a = a„. Note that for n — 9, we 
have a = 90 as one of the pitchfork bifurcation points. 
The Poincare sections displayed in Fig [2] shows that the 
stable channel orbit at a = 90 (Fig[2lja)) bifurcates and 
gives birth to two new stable orbits (Fig[2Ib)) while the 
channel orbit itself becomes unstable. Thus, pitchfork 

bifurcations take place at a„ = 2,12,30,56,90, and 

anti-pitchfork at an = 6,20,42,72, This can be ob- 
served in the plot of Tr J(a) as a function of a shown in 
Figl 



IV. QUARTIC OSCILLATOR STATES AND 
REDUCED DENSITY MATRIX 



A. Quartic oscillator spectra 

The quantum spectrum of the quartic oscillator is ex- 
tensively studied and reported [19, 22, 28]. For the pur- 
poses of this study, we note that two classes of eigen- 
states can be identified. The first one is what we call a 
generic state whose probability density |^'„(a;,j/)p cov- 
ers the entire accessible configuration space. Most of the 
eigenstates fall in this class and they are instances of 




FIG. 3: The linear stability of the channel orbit as a function 
of a. The orbit is stable for |Tr J(a)| < 2. The pitchfork 
bifurcation points are indicated by circles and anti-pitchfork 
bifurcations are indicated by squares. 



Berry's hypothesis that the Wigner function of a typical 
chaotic state condenses on the energy shell [1^. In Fig 
|31[a), we show the expansion coefficients for the 1973rd 
eigenstate of the quartic oscillator counted sequentially 
from the ground state for a = 90. Notice that the state 
is delocalizcd over a large set of basis states. These class 
of states are well described by random matrix theory. 
The second class of states is the localized states, which 
has enhanced probability density in the vicinity of the 
underlying classical periodic orbits. Theoretical support 
for this class of states based on semiclassical arguments 
is obtained from the works of Heller [2^ , Bogomolny and 
Berry [g^]. As a typical case, Fig[4l^b) shows the expan- 
sion coefficients for the 1972nd state which is localized 
over very few basis states in contrast to the one in Fig 
|4fa). In this work, we concentrate on a subset of such 
eigenstates whose probability density is concentrated in 
the vicinity of the channel periodic orbit. This set of 
states are nearly separable and can be approximately la- 
belled by a doublet of quantum numbers (N, 0) using the 
framework of adiabatic theory [H, [2^ . Note that such a 
labeling is not possible for the generic states since they 
are spread over a large number of basis states. 



B. Reduced density matrix 

In this section, we compute the eigenvalues of the RDM 
and the entanglement entropy of the quartic oscillator 
eigenstates as a function of the chaos parameter a. In 
terms of the expansion coefficients in Eq. ([7]), the ele- 
ments of RDM, Hx, can be written down as. 



M 

E 

ni — l 



K 



712 ^ni,n2 "'Til ; 



where the normahzation constant K.. 



1 if ni 



(9) 



n2 



and = 1/2 if rii ^ n2. 
has been traced out. 



with elements 



\n2\p 



{y)\ 



In this case, the y-subsystem 
Similarly another RDM, Rj,, 
-2) can be obtained by tracing 



over X variables. Let A represent the eigenvector ma- 
trix of order (Af -t- 2)/2 with elements a„^.„2, where 
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The exponential localization implies that a„ 



for 



5000 10000 

FIG. 4; Quartic oscillator eigenstates for a = 90 in the unper- 
turbed basis, (a) 1973rd state (delocalized) , (b) 1972nd state 
(localized). The inset is the magnification of the dominant 
peak. The eigenvalues of the RDMs for (c) 1973rd state and 
(d) 1972nd state. The inset in (d) is the magnification of the 
dominant eigenvalues that display degeracy. 



^1,712 — 0, 2,4....M labels the rows and columns respec- 
tively. Then, in matrix language, the RDM = A'^A 
is matrix of order (M -I- 2)/2. 

In our case, M = 318 and we numerically diagonal- 
ize the RDM of order 160. The eigenvalues of RDM for 
a typical delocalized state and a localized state is plot- 
ted in Fig|4l^c,d). In general, the dominant eigenvalues 
fall exponentially, though with different rates, for both 
the generic and typical localized state indicating that 
the Schmidt decomposition provides a compact represen- 
tation for the given eigenstate. Earlier such a behavior 
was noted for coupled standard maps 0. The first few 
dominant eigenvalues of RDM for localized states display 
(near-)degeneracy (see FigHl^d)). This arises as a conse- 
quence of (i) Civ symmetry of the potential due to which 
the eigenvector matrix is symmetric, i.e, a„i_„2 — a„2^„j 
and (a) the localization is exponential in the direction 
perpendicular to that in which the quanta of excitation 
is larger [28j . i.e, ajsi,n2 ^ exp(— [jn2), where > is a 
constant independent of N . 

The origin of near-degeneracy can be understood by by 
considering a simple model of 4 x 4 symmetric eigenvector 
matrix (the state number index n is suppressed such that 



/ 



ao,o 

02,0 
O-N.O 



02,2 

O.N,2 



O'N.2 

aN,N 



aN+2,0 
aN+2,2 
aN+2,N 



(10) 



aN+2,0 aw+2,2 aAr+2,JV aN^2,N+2 



Here we have only used the one-dimensional quartic os- 
cillator quantum numbers (0, 2, TV, iV -|- 2) because the lo- 
calized states can be approximately well represented by 
all possible doublets arising from these quantum num- 
bers. For instance, an adiabatic separation with the 
{N, 0) manifold gives a good estimate for the energy of 
its localized states The representation gets better 

as we add more ID quantum numbers to the list above. 



ni ~ n2- Further, a„j_„2 « if ni,n2 « N. Thus, 
elements aN,N ~ aN+2,N ^ 0.^+2^+2 ^ ao,o ~ 0. Then, 
we can identify a block matrix B with non-zero elements 
as, 

V o,N,2 aN+2,2 J 
Then, the eigenvector matrix P can be approximated as, 

Under the conditions assumed above, the RDM separates 
into two blocks which are transpose of one another. Thus, 
the reduced density matrix will have the form. 



R = P^P = 



BBT 
B^B 



(13) 



Since the eigenvalues remain invariant under transpo- 
sition of a matrix, i.e, the eigenvalues of BB"^ and 
B'^B are identical and hence we obtain the degeneracy. 
Though we use a 4 x 4 matrix to illustrate the idea, this 
near degeneracy would arise for any eigenvector matrix 
of even order, if the symmetry and exponential decay 
conditions are satisfied. 

For the localized state shown in Fig IH^b) , N — 264 
and the dominant eigenvalue of RDM using the approx- 
imate scheme in Eqns (|10j|13l) . is Ai = 0.4434. This is 
doubly degenerate and compares favorably with the ex- 
act numerical result of 0.4329. As observed in Figdl^d), 
the degeneracy breaks down as we travel down the index. 
As pointed out, the dominant eigenvalues of RDM corre- 
spond to definite ID quantum oscillator modes that ex- 
hibit exponential decay in the perpendicular mode. This 
is not true of all the oscillator modes and hence the de- 
generacy is broken. 



C. Entanglement entropy 

Entanglement entropy for each eigenstate is computed 
from the eigenvalues of the RDM using Eq In Fig 
[5l we show the entanglement entropy of the quartic os- 
cillator at a = 30 for one thousand eigenstates starting 
from the ground state. The localized states have values 
of entanglement entropy much lower than the local aver- 
age as seen from the dips in the figure. Most of them are 
much closer to zero and substantiate the fact that they 
are nearly separable states. In the next section we will 
show that the entanglement entropy of localized state is 
modulated by the bifurcation in the underlying channel 
periodic orbit. 

The generic delocalized states, on the other hand, form 
the background envelope seen in Fig O These chaotic 
states are not affected by the bifurcations in the iso- 
lated orbits. It is known that such delocalized states can 



5 





1000 



10 20 30 40 50 60 70 80 90 100 

a 



FIG. 5: (Color Online) Entanglement entropy for the quartic 
oscillator at a = 30 from ground state to 1000th state. The 
localized states have lower value of entanglement entropy as 
seen from the dips in the curve. The solid red curve is Srmt, 
the RMT average of entanglement entropy. 



FIG. 6: (Color Online) Entanglement entropy as a function 
of a. The three curves correspond to different (A*', 0) type 
localized states; solid circles (240,0), open circles (200,0) and 
squares (210,0). The positions of pitchfork bifurcations (tri- 
angle up) and anti-pitchfork bifurcations (triangle down) are 
marked on both the x-axes. 



be modeled using random matrix theory and hence the 
distribution of their eigenvectors follows Porter-Thomas 
distribution [2^ . The entanglement entropy can also be 
calculated based on RMT assumptions and it is known 
to be [3(|, Srmt — ln(7M) where 7 ~ l/\/e and M is 
the dimensionality of the reduced density matrix. In the 
case of quartic oscillator, the Hilbert space is infinite in 
dimension and we take M to be the effective dimension 
Meff of the RDM. One indicator of the effective dimen- 
sion of the state is the inverse participation ratio of the 
eigenstates. Based on this measure and due to Civ sym- 
metry of the quartic oscillator, we have M^^j: = D where 
D is the state number. Thus, the effective dimension of 
RDM is, Meff = VD. Finally, we get for the entangle- 
ment entropy, 

SRMT=H7Meff)^H-fVD). (14) 

In Fig[5l Srmt is shown as solid red curve and it correctly 
reproduces the envelope formed by the delocalized states 
while the localized states stand out as deviations from 
RMT based result, namely, Srmt- 

V. ENTANGLEMENT ENTROPY AND 
BIFURCATIONS 

In this section, we show the central result of the pa- 
per that the entanglement entropy is a minimum at the 
points at which the underlying periodic orbit undergoes 
a pitchfork bifurcation. As pointed out before, the lo- 
calized states of the quartic oscillator are characterized 
by the doublet {N, 0) and are influenced by the chan- 
nel periodic orbit. We choose a given localized state, say, 
with N = 200 and compute the entanglement of the same 
state, i.e, (200, 0) state as a function of a. The state that 
can be characterized by the doublet (200, 0) will be a lo- 
calized state at every value of a. The result is shown in 
Fig [5] as the curve plotted with open circles. The values 



of a at which the pitchfork and anti-pitchfork bifurca- 
tion takes is marked in both the horizontal axes of the 
figure as triangle-up and triangle-down respectively. For 
the purpose of easier visualization, they are connected 
by vertical lines. Notice that the entanglement entropy 
attains a local minima in the vicinity of every classical 
pitchfork bifurcation and it attains a local maxima near 
every anti-pitchfork bifurcation. As Fig [6] shows, similar 
result is obtained for two different localized states with 
{N = 210, 0) and {N = 240, 0). All these locahzed states 
are in the energy regime of highly excited states where 
the classical system is predominantly chaotic. The strik- 
ing similarity between the classical stability curve for a 
particular periodic orbit, namely the channel orbit, in 
Fig. ([3]) and the variation of the entanglement of the lo- 
calized state is to be noted. We have also numerically 
verified (not shown here) that a similar result is obtained 
in the case of another potential where pitchfork and anti- 
pitchfork bifurcations of the channel periodic orbit play 
an important role, namely, in V{x, y) = A-y^ + /3x^?/^, 
where (3 is the chaos parameter. 

At a pitchfork bifurcation, as shown in Fig. [2l the 
fixed point corresponding to the channel periodic orbit 
loses stability and becomes a hyperbolic point. The cen- 
tral elliptic island seen in Fig[5fa), breaks up into two 
islands. The localized state that mainly derives its sup- 
port from the classical structures surrounding the stable 
fixed point suffers some amount of delocalization, but 
is largely supported by the stable regions. At an anti- 
pitchfork bifurcation, the hyperbolic point becomes an 
elliptic fixed point and the orbit has gained stability and 
a small elliptic island just comes into existence. Hence, 
the eigenstate is still largely delocalized since the small 
elliptic island is insufficient to support it. This heuristic 
picture which is quite sufficient to explain oscillations in 
localization measures is seen to be surprisingly valid even 
for the somewhat less intuitive measure of entanglement. 
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It has been noted earlier that when the corresponding 
classical system undergoes a pitchfork bifurcation, the 
entanglement entropy defined by Eq (|4]) attains a maxi- 
mum [9|, and it has been conjectured to be a generic prop- 
erty. We note that this, apparently contradictory result, 
is however in the context of an equilibrium point under- 
going a bifurcation and the relevant state is the ground 
state, whereas in the case we are studying here the or- 
bit that is bifurcating is a periodic orbit and the states 
are all highly excited. In this situation there is a much 
tighter correlation between more conventional measures 
of localization (Shannon entropy, participation ratio etc.) 
and entanglement. 

As the parameter a is increased, the quartic oscillator 
gets to be predominantly chaotic and this should imply 
increase in entanglement. However, this is true only for 
the generic delocalized states as seen in Fig [5] The local- 
ized states are influenced not so much by the increasing 
volume of chaotic sea but by the specific periodic orbits 
that underlie them. Hence, for these states, it is only 
to be expected that the qualitative changes in the phase 
space in the vicinity of the corresponding periodic orbits 
affect quantum eigenstate and hence its entanglement as 
well. This can be expected to be a generic feature of 
entanglement in quantum eigenstates of mixed systems. 



VI. CONCLUSIONS 



as a bipartite system. We study the effect of classical 
bifurcations on the entanglement of its quantum eigen- 
states. The quartic oscillator is a classically chaotic sys- 
tem. One particular class of eigenstates of the quartic 
oscillator, the localized states are scarred by the chan- 
nel periodic orbits. We have shown that the entangle- 
ment entropy of these localized states is modulated by 
the bifurcations in the underlying channel periodic orbit. 
When this orbit undergoes a pitchfork bifurcation, the 
entanglement attains a local minimum and iwhen it un- 
dergoes an anti-pitchfork bifurcation the entanglement is 
a local maximum. Physically, this is related to the pres- 
ence or the absence of elliptic islands in the phase space 
in the vicinity of the channel orbit. We expect this to 
be a general feature of bipartite quantum systems whose 
classical analogue display bifurcation features. 
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